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ABSTRACT 

We give an introduction to the problems faced on the way to a rehable lattice QCD 
computation of B-physics matrix elements. In particular various approaches for 
dealing with the large scale introduced by the heaviness of the b-quark are mentioned 
and promising recent achievements are described. We present perspectives for future 
developments. 
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1 B-physics and lattice QCD 

The truly beautiful results from recent B-physics experiments [TJ |2l El IH E] repre- 
sented highlights of this conference. Some of them, such as the B — B mass difference 
Am-d, require knowledge of QCD-matrix elements for their interpretation in terms 
of parameters of the standard model of particle physics and its possible extensions. 
This motivates investigations in lattice QCD, our best founded theoretical formu- 
lation of QCD. Indeed, this formulation allows for the computation of low energy 
hadronic properties through the Monte Carlo evaluation of the Euclidean path inte- 
gral. While such a computation necessarily involves approximations, which we will 
discuss, the important property of the lattice approach is that all approximations 
can be systematically improved. 

Before going into the details we summarize the goals of lattice QCD com- 
putations with b-quarks. They motivate the considerable effort involved. 

• The determination of parameters of the CKM matrix, which in the Standard 
Model are fundamental parameters of Nature. In particular the unitarity tri- 
angle should be determined and over-constrained in order to test the Standard 
Model. The necessary non-perturbative matrix elements should be computed 
through lattice QCD. 

• A precise computation of the b-quark mass, which enters many phenomeno- 
logical predictions and plays a role in grand unification and other questions 
beyond the Standard Model. 

• A determination of the spectrum and lifetimes of b-Hadrons. 

• Non-perturbative tests of the Heavy Quark Effective Theory (HQET) which 
is applied frequently to simplify the dynamics of heavy quarks, but is difficult 
to test experimentally. 

The starting point of a lattice computation is the QCD Lagrangian, formu- 
lated on the discretized 4-dimensional Euclidean space-time, i.e. on a hyper-cubic 
lattice with spacing a [HI El IE] • The beauty of this theory is that it contains only 
the (bare) gauge coupling go and the (bare) quark masses as parameters. After 
these parameters have been fixed by a small set of experimental observables, say the 
masses of proton, pion, kaon, D-meson and B-meson, all observables become pre- 
dictions of the theory^. In order to have a finite number of variables in the Monte 
Carlo evaluation of the path integral, one considers a finite space time of linear size 
L, mostly with periodic boundary conditions. One then has to approach 

^For simplicity we neglect the top quark and assume that electroweak effects are treated as 
perturbations. 
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=^ the infinite volume limit, L ^ oo and 

=^ the continuum limit, a — > . 
While the statistical errors of the Monte Carlo decrease oc computer time at 
fixed a, L, a reduction of the total error requires to approach the above limits. There- 
fore the total error decreases much more slowly as computers get faster. Progress is 
made by better formulations of the problems, the development of better computa- 
tional algorithms and the rapid increase in computer speed. Here we want to explain 
the particular challenges one faces in reaching small overall errors in B-Physics and 
discuss recent advances in facing them. Some results will be shown, mainly to illus- 
trate the progress that has been made. A more comprehensive list of results can be 
found in [giin]. 

2 The challenges 

2.1 Renormalization 

One of the challenges that one faces has to do with the fact that interesting transi- 
tions originate from the electroweak interactions, which we can not treat simultane- 
ously with QCD in the simulations. One then adopts the (good) approximation to 
treat the electroweak interactions at the lowest non-trivial order in the electromag- 
netic and weak coupling. Consequently one has an effective Hamiltonian, valid at 
energies far below the electroweak scale, which contains the quark and gluon fields 
but not the photon and the electroweak bosons. For example a left-left 4-fermion 
operator is one of the operators in the effective Hamiltonian and its matrix element 
determines Arrid- The renormalization of such operators is a non-trivial task, but 
besides perturbation theory [TT] powerful non-perturbative approaches have been 
developed [121 CHI Ull HHI and they are continuously improved and applied to new 
operators fHl- We will come back to this challenge in our discussion of the HQET. 

2.2 The multi-scale problem 

Even after eliminating the electroweak scale by considering low energy matrix ele- 
ments, b-physics always contains two more scales apart from the typical QCD scale 
A = O(lGeV). The first is the scale of the light quark masses, or better the mass 
of the pion and the second is the large mass of the b-quark itself. To understand 
what this means in practical terms, consider a lattice with 32^ x 64 points, as it is 
presently possible in the quenched approximation but still prohibitively expensive 
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for the full theory.^ Choosing L 2.5 fm and a ^ 0.1 fm, one notices that L is 
comparable to the Compton wave length m~^, namely L ^ and at the same 

time the b-quark mass is larger than the inverse lattice spacing: mb ~ Instead 
we should have 

L » , (1) 
a < m^^ , (2) 

in order to keep finite size effects due to pion propagation around the periodic world 
small ( eq. (P)) and to properly resolve the propagation of a b-quark (eq. ©)• It 
is known that effects of the first type will become rapidly (exponentially) small 
when L > 4m~^ ITH] . while the dominant discretization errors due to a heavy 
b-quark are given by 0((amb)^) if the theory is 0(a)-improved (i.e. linear a effects 
are removed [121 1201 )• The latter behavior sets in roughly for a < \m^^ fH] . 
On a 32^ x 64 lattice one is far from satisfying these constraints. An additional 
very relevant factor is that the computational cost of simulating full QCD grows 
rapidly when quark masses are lowered. This means that presently simulations at 
the physical values of the light quark masses are impossible. 

In short: the light quarks are too light and the heavy quarks are too heavy 
for present computing capabilities. 

In order to obtain results, a reformulation of the theory (or specific prob- 
lem) is needed or extrapolations have to be performed. 

2.3 Chiral extrapolations 

Reliable extrapolations of numerical data are possible when sufficient analytic knowl- 
edge of the functional dependence on the extrapolation parameter exists and when 
the data are available in a range where the analytic formulae apply. Concerning the 
light quark mass dependence, chiral perturbation theory furnishes an expansion in 
terms of them, which is also applicable in the case of heavy-light mesons [221 
Unfortunately it remains unclear up to now, how small the quark mass has to be in 
order that this expansion is applicable at the quantitative level. Most current simu- 
lations reach light quark masses which are about half as heavy as the strange quark 
mass and an agreement with the analytic formulae could not yet be established |24j . 

^After integration over the quark fields the Boltzmann weight of the path integral contains a 
factor of the determinant of the Dirac operator. In the quenched approximation the dependence 
of this determinant on the gluon fields is neglected. In a perturbative language this corresponds 
to neglecting quark loops, but keeping all the gluon exchanges. 
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For recent reviews of the numerical situation with an emphasis on B-physics see 
refs. (23 123- 

As an example we would like to mention only that estimates of the un- 
certainty for the phenomenologically important ratio ^ = (Fb^ V^bJ / {Fb^ \/^b^) 
range approximately from 5 % [2Z| to 10 % [SH]. In our opinion these estimates have 
to be substantiated by simulations with smaller quark masses, where contact with 
the chiral perturbation theory expansion can be demonstrated. It appears that this 
requires both faster computers and the development of algorithms which perform 
better, in particular at small quark masses. First steps have been taken [^1^011^ . 
and new ideas exist |22|- These very important developments are not specific to 
B-physics. We therefore do not discuss them further. 



2.4 Heavy quarks on the lattice. 



Continuum extrapolation 



0.04 



Before coming to the methods for 
dealing with the problem of a heavy 
b-quark, let us give an illustration of 
the extent of the problem. It has 
been very well investigated for the 
charm quark E31 , which is a fac- 
tor 4 lighter than the b-quark. In 
Fig. we show the renormalization 
group invariant (RGI) charm mass^, 
computed on four lattices, between 
a = 0.1 fm and a = 0.05 fm. Three 
definitions of the renormalized quark 
mass were considered, differing by dis- 
cretization errors O(a^). The figure 
shows that on lattices of the chosen force between static quarks. 
size (from 16^ x 32 to 32^ x 64) one 

can control the discretization errors by extrapolation, but it is also obvious that this 
will not be possible for much heavier quarks. For b-quarks other approaches have 
to be considered. The following list summarizes what is being discussed. 

1. Extrapolation in the mass of the heavy quark, rrih- 

2. Effective theories which implement an expansion in 1/m-t,: NRQCD|35j and 
HQETlSn]. 




Figure 1: RGI charm mass from 133]. The 
length scale tq 0.5 fm is derived from the 



definition is given in eq. below. 
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3. The "Fermilab approach" ^7\ . 

4. Anisotropic lattices with atimc c^space |HH1 1^^] • 

5. Combinations, in particular of 1. and 2. 

6. Extrapolation of finite volume effects in the quark mass I41j. 

1. is clearly applicable once rrih can be made large enough for the functional form, 
a power series in l/rrih, to be accurate. If F{mh,amh) is the desired observable at 
finite lattice spacing and (unphysical) heavy quark mass, one has to evaluate 

F(mb, 0) = lim lim F(mh, arrih) (3) 

by two subsequent extrapolations whose order is important. In practice a residual 
uncertainty remains due to the assumptions made in the extrapolation rrih — >■ mb. 

2. Effective theories should be rather accurate since A/mb is a small expansion 
parameter, which controls both NRQCD and HQET. For reasons which will become 
clear below, NRQCD has predominantly been used in recent years although its 
continuum limit does not exist. In this theory one must keep the lattice spacing 
finite and control the discretization errors by adding terms to the Lagrangian that 
remove them approximately. An additional source of errors is that the coefficients 
in the Lagrangian are usually determined perturbatively. Below we will report on 
recent progress in HQET on the lattice. 

3. The Fermilab approach was discussed at last year's conference |26j . 

4. The possibility whether anisotropic lattices do permit a situation aspace^^b > 1 
for the interesting matrix elements is still under discussion In any case, 
giving up the symmetry between space and time allows for dimension four operators 
in the Lagrangian which break Euclidean invariance. To obtain a Euclidean invariant 
continuum limit (and thus Lorenz invariance after rotation to Minkowski space), 
these parameters have to be tuned properly. For full QCD with dynamical fermions, 
this is a non-trivial task if non-perturbative precision is desired. 

5. Combining the effective theory at the lowest order, (the "static theory"), with (1.) 
turns extrapolations into an interpolations, reducing the influence of assumptions 
concerning the l/m/j-dependence considerably. 

6. This new idea put forward at Tor-Vergata j^E] will be explained below. 

The dominant approach in the last decade has been to compare several 
methods with their strengths and weaknesses and apply the result to phenomenology 
if different methods agree. New developments offer the chance to establish precise 
results without recourse to crosschecks through other methods. 
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3 New developments 

3.1 Non-perturbative HQET 

The lattice Lagrangian of (velocity zero) HQET, 

>Chqet = V'h Do A - I^V'hD^^/'h - -^^hB ■ crtph + ■■■ (4) 

has the same form as the continuum one to order 1/mb; only the definitions of 
the covariant derivatives and the chromomagnetic field strength B are of course 
lattice specific. Together with a similar expansion of the operators who's matrix 
elements one is interested in, it implements a systematic expansion in terms of 
1/mb for B-mesons at rest [32] . Despite its attractive feature of having a continuum 
limit order by order in the 1/mb-expansion, it has not been applied very much in 
recent years. The reason is threefold. 

First, already in lowest order of the effective theory, called the static ap- 
proximation, statistical errors grow rapidly as the Euclidean time-separation of cor- 
relation functions is made large (Fig. |21 filled symbols). But it is in the large time 
range, say a;o > 1.5 fm, where masses and low energy matrix elements may safely be 
extracted. 

Second, the number of parameters in the effective theory grows with the 
order in the expansion in 1/mb. 

Third, these parameters have to be determined non-perturbatively; other- 
wise the continuum limit does not exist |33] (fine tuning of parameters). This fact is 
due to the mixing of higher dimensional operators, such as ^^ip^T^'^ip^^ with lower 
dimensional ones, such as ip-^Doipii. A perturbative estimate of the parameters in 
the Lagrangian (or equivalently of the mixing coefficients) to order Qq'' would leave 
a perturbative 

remainder ~ ^0^'+^^ ~ a'^ [ln(aA)]-('+^) ^ oo . (5) 

This basic fact is unavoidable in an effective theory formulated with a cutoff. 

Recently it was shown that the first point is overcome by considering al- 
ternative discretizations of the static theory ^3], which differ only in the way the 
gauge fields enter the latticized covariant derivative. So-called HYP-links jlSJ US] 
correspond to the triangles in Fig.|21and result in errors at the % level at xq = 1.5 fm 
with an only slow growth cLS Xq is increased. Very importantly, it was also shown 
that a-effects with this new discretization are small OI37|. 



7 



The second and third point 
above can be solved in one go if 
the parameters of HQET are non- 
perturbatively determined from those 
of QCD. In this way the predictive 
power of QCD is transfered to HQET. 

The basic idea how to do 
this 03 CEl, illustrated in Fig. El is 
easily explained. In a finite volume 
of linear extent Lq = 0(0.2 fm), one 
may realize lattice resolutions such 
that -C 1 and the b-quark can 

be treated standard relativistic 
fermion. At the same time the energy 
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Figure 2: Noise to signal ratio of a B-meson 
correlation function in static approximation 
at a ^ 0.08 fm JJ^L The red bullets are for 
the original Eichten-Hill action while purple 
triangles correspond to the alternative dis- 
cretization using HYP-links. 



scale 1/Lq = 0(1 GeV) is still signifi- 
cantly below TTib such that HQET applies quantitatively. Computing the same suit- 
able observables in both theories relates the parameters of the HQET Lagrangian 
to those of QCD (matching). One then remains in HQET, changing iteratively to 
larger and larger volumes, and computes HQET observables in each step. Finally, 
in a physically large volume (linear extent 0(2 fm)) the desired matrix elements 
are accessible. At the end any dependence on the unphysical intermediate volume 
physics is gone except for terms of the order 0((Lomb)~*-"''''^-*) if the effective theory 
was considered up to order n. 

The strategy is formulated in such a way that the continuum limit can be 
taken in each individual step. To explain this we take a look at the simple equation 
which - at the lowest non-trivial order in l/rrib (static approximation) - relates the 
B-meson mass to the mass of the b-quark. For definiteness we take the RGI quark 



experiment 



Lattice with arrib <^ 1 



me = 5.4 GeV 



{L2 



■ stat 



(Tm('Ui) 



Li — 2*Lo 

rstat(-^^l) 



(Tm('Uo) 



r(Lo,M) 



stat 



Figure 3: Connecting experimental observables to renormalized HQET. Fgtat is a 
renormalized quantity in HQET and a^{g'^{L)) connects rstat(-^) o-nd rstat(2-L). In 
the chosen example, the experimental observable is the mass of the B meson. 
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mass. It is given by the large /i asymptotics of the running mass, m(/i) (in any 
scheme) via, 

Mb = Jim mb(/^) [26o^'(/i)J , (6) 

with h^^df) the lowest order coefficients of the beta-function and the anomalous 
dimension of the quark mass, respectively (conventions as in ^15,). In contrast to 
mb(/i), the RGI-mass, Mb, is scheme independent. 

In static approximation it is related to the mass of the B- meson, me, via 

me = fitat - rstat(iv2) + rstat(i^2) " r,tat(^o) + r(Lo, Mb) +0(AVMb) . (7) 
a^o in HQET a-^o in HQET a-*o in QCD 

Here r(L,Mb) denotes the energy of a state with quantum numbers of a B-meson 
but defined in a finite volume world of linear extent L. The exact definition of 
this state |1H] is not important to understand the idea but is quite relevant for 
the success of a numerical computation of Mb. rstat(-^^) is the same energy but 
evaluated in static approximation and i^stat denotes the energy (mass) of a B-meson 
state in large volume in static approximation. As mentioned before we have Lq ~ 
0.2fm, L2 = 4Lo. 

To appreciate eq. one should first note that energies in the effective 
theory are related to energies in QCD by a shift mbare which is universal in the 
sense of being independent of the state. This corresponds to a term mbare iph '4'h 
in eq. (j3)), which we have dropped there following standard conventions. Since the 
operator ip^i Dq ipt^ mixes with the lower dimensional one ip-^ tp^ under renormaliza- 
tion, the parameter mbare is linearly divergent (~ 1/a) and must be determined 
non-perturbatively. Its universality means 

TUB = -Estat + "Ibare , (8) 

r(L,Mb) = rstat(^) + ^bare ; (9) 

with one and the same mbare at fixed a. One may now use eq. to determine 
the parameter mbare in the effective Lagrangian from QCD and then insert it into 
eq. (jHl) to determine me. This represents the general logics for obtaining results in 
the HQET. In order to arrive at the continuum limit of the prediction, one groups 
terms as in eq. ((Tj), where mbare drops out of the energy differences i^stat — rstat(-Z^2) 
and rstat(-^^2) — rstat(-^o) and the continuum limit can be taken separately for each of 
the terms as indicated. Indeed, in the quenched approximation a precise continuum 
limit has been taken for all terms except for the energy difference ii^stat— rstat(-^^2)- 
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Figure 4: The combination AE = E'stat — rstat(-^2) as a function of the lattice spacing. 
Results with the new discretizations '^441 (r.h.s) are compared with old results \4^ - 

Here flH] uses values for i^stat from the literature (l.h.s. Fig. 0)). With the new 
discretization [H] more precise results are obtained (r.h.s. Fig. H}, and they also 
have the linear 0(a)-errors removed. The full analysis with the new data has not yet 
been performed. Using for the time being the value of Lo(£'stat — rstat(-^2)) = 0.46(5) 
one obtains jUj (we do not distinguish 0(A/(LoMb)) and O(AVMb) ) 



roMb = 16.12(28) + O(AVMb) 



_MS/;^MS 



m 



4.12(8)GeV + 0(AVMb) (10) 



in the quenched approximation and with tq = 0.5 fm. It is worth emphasizing that 
such a result is based on the non-trivial relation between the bare quark masses on 
the lattice and the RG I- masses established in [T^ I50j. 

3.2 Results for Fb^ 

As a further new development, we show in Fig. |31a recent computation of the decay 
constant of the Bs meson in quenched approximation using method number 5 
in our list of Sect. 12.41 HQET predicts the mass dependence 



i^pSv^/Cps(M/A) = <l>RGi + 0(l/mps) 



of the decay constant Fps, where $rgi is independent of the heavy quark mass. It 
has recently been computed in the continuum limit of the static approximation |44j . 
The factor Cps is a function of the ratio of the RGI mass of the heavy quark and 
the QCD Lambda parameter. It is now known quite accurately from perturbation 
theory due to the 3-loop result of Also the numbers of Fps at finite quark 

mass shown in Fig. El have been extrapolated to the continuum limit jHHEll- While 
the subsequent interpolation in 1/mps is obviously safe, an extrapolation with just 
results for l/(romps) > 0.15 would depend on the functional form assumed. The 
present result F^^ = 206(10) MeV [HI] is in good agreement with the recent one from 
the Tor-Vergata group 
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Figure 5: The continuum limit quenched pseudo scalar heavy- strange decay constant 
as a function of the inverse heavy-strange meson mass, mps l.^ij /. 

3.3 The Tor-Vergata approach: 

Extrapolation of finite volume effects in the quark mass 

The starting point of this method is the same as in Sect. 13.11 in an intermediate 
volume (e.g. of size Lf x 2Li) the lattice spacing can be made small enough to be 
able to treat b-quarks as relativistic fermions. The decay constant F{Mh, Li) can 
then be computed for = Mb. While in the non-perturbative HQET, this was 
used to obtain the parameters in the effective theory, de Divitiis et al. remain in the 
relativistic theory, but introduce as their central observables finite size effects ap in 
the following way'* 



The idea is that the finite size effects crp depend strongly on the dynamics of the 
light quark, but if Mh is large enough, they hardly depend on that variable. It is 
thus expected that finite volume effects can smoothly be extrapolated in the heavy 
quark mass, dp (Mb, Li) = limM^^Mt (TpiMh, Li). 

Setting Li ^ 0.4 fm, a numerical computation was performed in the quenched 
approximation. Of course, the three factors in eq. (fT^ are obtained by an extrap- 
olation to the continuum limit at fixed M/i, followed by an extrapolation in 
to Mb for (Tf, while F(Mb,Li) could be determined directly at Mh = Mb. The 
most difficult extrapolations are those concerning the step 2Li — > 4Li, since there 

^For 4Li w 1.6 fm further finite size effects are expected to be very small (at least in the 
quenched approximation) . 



F(Mb,oo) ^F(Mb,4Li) 



F(Mb,4Li) ^ F(Mb,2Li) 

F(Mb,2Li) F{M^,L,) 
' ' 



xF(Mb,Li). 



(12) 



=o-F(Afb,2Li) =crF(A/b,Li) 
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(a/rJ2 1/M, 

Figure 6: Quenched computation of ap{M]^,2Li) with the light quark mass equal to 
the strange quark mass. 

a > 4Li/32 ^ 0.05 fm on lattices with maximally 32'^ x 64 points. They are shown 
in Fig. El 

Comparing the r.h.s. of Fig. IHl with Fig. one notices that (looking at the 
central values) in both cases the last point of the finite mass numbers differs from 
the result at mps = me^ by about 7%. 

The final result is |lTj 

Fb, = 192(6)(4)MeV. (13) 

Applying the same method also to the computation of the b-quark mass, yields [iU] 
m^(m^) =4.33(10) GeV. 

3.4 Summary 

The new developments discussed in this section attack the problem of a heavy quark 
mass in ways where all necessary renormalizations are performed non-perturbatively 
(only in Cps an 0{a{m-b)^) uncertainty remains) and the continuum limit can be 
taken in each step. They have been shown to be applicable in numerical compu- 
tations, yielding good accuracy after propagating errors through all the steps. In 
fact, all the errors quoted in this section can be further reduced in the quenched 
approximat ion . 

4 Perspectives 

We start with some possible straight forward improvements of the above results. 
First, the precision in Sect. 13.11 and 13.21 can easily be improved. Once this is done, 
one may quantify the size of corrections rather precisely. In this context we 

note that it is surprising that so far there is no sign of (1/mc)^ terms in the charm 
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region (see Fig. Eland Fig. Better precision is desirable to pin those down. 
Second, a very promising improvement would be to combine the method of Sect. 13.31 
with Sect. 13.11 the finite size effects ap can be computed in the HQET. Then the 
mass-extrapolation in Fig.lHlwill also be turned into an interpolation and even larger 
precision and confidence can be achieved. 

In our opinion it is most important (and more ambitious) to extend the 
calculations in the HQET to the level of including 1 / m-b corrections. We do not see 
any major obstacles on the way to reaching this goal. The severe problem of power 
divergences can be overcome by the non-perturbative matching of QCD and HQET 
sketched in Sect. 13.11 In contrast to the other methods discussed, the consequent use 
of an effective theory to describe the b-quark is a way to entirely eliminate the mass 
scale of the b-quark and thus the necessity of using large lattices. This then opens 
the possibility to obtain results in full QCD, once the usual (and large) problems 
with dynamical fermions, which have nothing to do with the heaviness of the b, 
are solved. Indeed, as emphasized in Sect. 12.31 we finally do need results with light 
dynamical quarks in order to apply them to phenomenology with confidence. 
Acknowledgments. I am indebted to my collaborators in the ALPHA-coUaboration 
for the possibility to present new results at this conference, prior to publication. I 
thank the Tor Vergata group for communications on their approach and for sending 
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